Infrared spectroscopy of electronic bands in bilayer graphene 
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We present infrared spectra (0.1-1 eV) of electrostatically gated bilayer graphene as a function of 
doping and compare it with tight binding calculations. All major spectral features corresponding 
to the expected interband transitions are identified in the spectra: a strong peak due to transitions 
between parallel split-off' bands and two onset-like features due to transitions between valence and 
conduction bands. A strong gate voltage dependence of these structures and a significant electron- 
hole asymmetry is observed that we use to extract several band parameters. Surprisingly, the 
structures related to the gate-induced bandgap are much less pronounced in the experiment than 
predicted by the tight binding model. 



Since the first successful attempt to isolate graphene 
[l[, this two-dimensional material remains in the focus 
of active research motivated by a unique combination of 
electronic properties and a promising potential for appli- 
cations Q. Its infrared response, like many other trans- 
port and spectral properties, is notably distinct from the 
one of conventional metals and semiconductors. For ex- 
ample, the optical conductance Re G{u}) of monolayer 
graphene, which describes the photon absorption by a 
continuum of electronic transitions between the hole and 
electron conical bands, remains constant in a broad range 
of photon energies and equal to Gq = e'^/Ah 0, 0, Si- 
Quite remarkably, the optical transmittance of single car- 
bon layer in this range depends solely on the fine struc- 
ture constant In bilayer graphene, where the inter- 
layer electron hopping results in two extra electron and 
hole bands separated from the main bands by about 0.4 
eV, one expects to see a set of intense and strongly dop- 
ing dependent infrared structures 0, 0, 0] sensitive to 
various band details and quasiparticle scattering rates. 
This makes infrared spectroscopy a powerful probe of the 
low-energy electronic dispersion in graphene, especially 
in combination with apossibility to electrostatically con- 
trol the doping level [s, E, El ■ Here we present infrared 
spectra of bilayer graphene crystals in a broad doping 
range, which allows us to observe new features, in partic- 
ular a significant electron-hole asymmetry. By comparing 
data with the tight binding Slonczewski-Weiss-McClure 
(SWMcC) model p^] we identify interband transitions 
and determine some band parameters. 

Bilayer graphene is considered to be particularly 
important for electronics applications by virtue of a 
bandgap that opens when there is difference between the 
electrostatic potential of the two layers [l^ [3, [l^ [l^, 
17 , [Tl |. Angle-resolved photocmission (ARPES) mea- 



surements indicate such a gap in potassium doped bilayer 
graphene epitaxial ly g rown on SiC [3] ■ Although trans- 
port experiments 



17, m 



demonstrate that a bandgap 
also opens in gate-tunable bilayer graphene flakes, no 
spectroscopic information about the size of the gate- 
induced gap is currently available. The analysis of in- 
frared data allows us to get further insight into this issue. 
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FIG. 1: (a) Schematic view and a micrograph of the used 
bilayer graphene device, (b) Infrared reflectance of graphene 
flake (blue solid line) and of bare substrate (red dotted line) 
(taken at T = 10 K and Vg = -1-100 V). Left inset: Bernal 
stacking of bilayer graphene and relevant hopping terms, right 
inset: resistivity at 10 K as a function of the gate voltage. 
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The sample used in this study is a large (~100 /im) bi- 
layer graphene flake (Graphene Industries Ltd.) on top 
of an n-doped Si substrate covered with a 300 nm layer 
of Si02 (Figlll^a)). A field-effect device configuration 
allowed us to simultaneously measure the DC resistiv- 
ity and infrared reflectance as a function of the applied 
gate voltage Vg. Optical spectra in the photon energy 
range 0.1-1 eV were collected at « 10 K with an infrared 
microscope (Brukcr Hyperion 2000) focussing the beam 
on a spot of about 30 microns in diameter. The ab- 
solute reflectance of graphene, i?fiakc, and of the bare 
substrate, i?oxidc, (FiglUb)) were obtained by using a 
circle of gold deposited close to the sample as a refer- 
ence mirror. The bare substrate spectrum features in- 
tense optical phonon modes in Si02 below 0.15 eV and 
a dip at 0.7 eV due to the Fabry-Perot effect in the Si02 
layer. The change of the absolute reflectivity introduced 
by graphene AR — Ruakc — -Roxide is small but repro- 
ducibly measurable as we checked on a second sample. 
By taking difference spectra, we largely cancel spurious 
optical effects such as a weak 0.4 eV absorption band due 
to some frozen water. The resistivity maximum that cor- 
responds to zero doping (FiglTJb), inset) is found to be 
at Vgo = -25 V instead of V, which we attribute to a 
charging effect by contaminant molecules. 

The curves of Ai?(w) between 0.2 and 0.6 eV are shown 
in Fig[2ja) as a function of the gate voltage from -100 V 
to -1-100 V. The spectra in this region are very sensitive 
to the gate voltage and show a significant asymmetry be- 
tween the electron (V^ > V^o) and the hole {Vg < Vgo) 
doping. Since it is more convenient to discuss the data in 
terms of the real part of the optical conductance G{uj), 
we extracted this quantity by a Kramers-Kronig (KK) 
constrained inversion [l^ of the raw reflectivity data. 
Due to a sensitivity of the inversion procedure to the 
systematic uncertainty (~ 0.005) of AR and to the data 
extrapolations beyond the experimental spectral range 
(we used graphite optical data Q as the most reasonable 
extrapolation) the inverted function Re G{uj) is likely to 
contain a spectrally smooth background as compared to 
Re G(w). Although this background does not allow us 
to determine accurately the absolute conductance, it af- 
fects the positions of spectral structures and their doping 
dependence to a much lesser extent. 

The spectra of Re G{u!) (Figlll^b)) reveal a prominent 
peak centered between 0.35 and 0.4 eV, whose intensity 
increases with the absolute value of the gate voltage and 
vanishes as Vg approaches V^o- Based on previous theo- 
retical works 0,M @] as well as on the calculations de- 
scribed below we assign this peak to a transition between 
the hole bands 1 and 2 (marked as C in FiglHl^e)) for 
Vg < Vgo and to the one between the electron bands 3 
and 4 (marked as B) for Vg > Vgo. The doping induced 
shift of the Fermi level away from the Dirac point ex- 
pands the momentum space, where this transition is al- 
lowed by the electronic occupation of the initial and the 
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FIG. 2: (a) Mid-infrared spectra of AR at T ?a 10 K as a 
function of the gate voltage Vg. The curves are separated by 
0.005; the dashed line is the zero level for the -1-100 V curve, 
(b) Real part of the infrared sheet conductance of bilayer 
graphene G{uj), derived from the reflectance curves (Fig. 2a) 
using a Kramers-Kronig inversion. The curves are separated 
by 0.5Go. Note that G{uj) possibly differs from the true con- 
ductance G{uj) by a spectrally featureless gate-independent 
background, as explained in the text. The dashed line is the 
correction (shown relative to the -1-100 V spectrum) used to 
generate Fig. |3]3. 

final states, and therefore increases the infrared intensity 
of the peak. 

The energy of this peak is given by the band separa- 
tion and is close to the interlayer vertical hopping param- 
eter 7i (shown in the inset of Fig[ljb)). In the case of 
precisely symmetric electron and hole bands, one would 
expect the same peak position for the positive and neg- 
ative gate voltages. However, the data reveal a clear 
asymmetry: at positive voltages the maximum (marked 
with red circles in Fig[2fb)) is higher in energy and shows 
a much stronger dependence on Vg than at negative volt- 
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ages (blue circles). As was pointed out in Refl22l the en- 
ergy of the peak on the electron and hole side taken close 
to the charge neutral point (Vgo = -25 V in our case) 
is equal to 71 + A and 71 — A respectively, where the 
parameter A is the potential difference between carbon 
sites A and B. These values in our case are 0.393 ± 0.005 
eV and 0.363 ± 0.005, which yields 71 = 0.378 ± 0.005 
eV and A = 0.015 ± 0.005 eV. The value of 71 is very 
close to 0.377 eV found in graphite However, it 



is somewhat smaller than 0.404 eV reported in Refl22 
for bilayer graphene flake. This suggests that the inter- 
layer distance, to which 71 is the most sensitive, may 
change from sample to sample. As far as A is concerned, 
there is much less agreement on the value of this pa- 
rameter in graphite in the literature. While the mag- 
netoreflection and de Haas - van Alphen measurements 
suggest that A is -0.008 eV (see Refill and references 
therein), infrared data [1, [2lj give a value of -1-0.04 eV. 
Our value agrees in sign with the infrared-based estimate 
in graphite, but is about 2-3 times smaller. This differ- 
ence can be understood using electrostatics arguments. 
In Bernal stacked graphite, each carbon layer is sym- 
metrically surrounded by two other layers, in contrast to 
bilayer graphene. Therefore one may expect the differ- 
ence between the (screened) Coulomb potential on sites 
A and B induced by charges on other layers to be larger 
in graphite. 

In order to get further insight, we compare the exper- 
imental data with calculations based on the tight bind- 
ing SWMcC model that proved to be very successful in 
graphite @, [13, The hopping terms considered are 

shown in the inset of Figlljb). The values of all band pa- 
rameters except 71 and A, which were determined above 
were taken from Ref. S 70 = 3.12 eV, 71 = 0.378 eV , 
73 = 0.29 eV, 74 = 0.12 eV and A = 0.015 eV. Note that 
they agree well with the values determined in Ref. [2^1 
using Raman spectroscopy. The doped charge and the 
Fermi energy can be directly determined for any given 
gate voltage using the known capacitance of the Si02 
layer [l3|. The Kubo formula was used to calculate G{uj) 
that was eventually Gaussian-broadened by 0.02 eV, in 
order to match the observed line widths. The reflectivity 
spectra were computed based on Fresnel equations using 
the known optical properties of the Si02/Si substrate. 
We begin with a calculation which assumes that the only 
effect of applying gate voltage is to shift the chemical po- 
tential and does not include the gate-induced bandgap. 

In panels (a) and (c) of Fig[31 the color plots of ex- 
perimental and calculated spectra of AR{uj, Vg) are rep- 
resented. One can notice a quite good correspondence 
between the energy and the gate voltage dependence of 
the strong spectral features. Having found that such 
an agreement is present in the raw reflectivity data, we 
proceed with a detailed experiment-theory comparison 
in terms of the optical conductance (FiglSlJb) and (d)). 
In view of the mentioned possibility that the extracted 



conductance curves contain a spectrally featureless back- 
ground, here we subtract from all spectra the same, 
i.e. gate- voltage independent, smooth curve shown as 
a dashed line in Figl^Jb). This curve is chosen in such 
a way that the corrected Re G(w, Vg = 100 V) coincides 
with the theoretical values in the regions around 0.2 eV 
and 0.6 eV, where no sharp structures are expected. 
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FIG. 3: (a) and (b): color plots of the raw AR{uj) and the de- 
rived Re G{ui) spectra as a function of and Vg. (c) and (d): 
AR and Re G{lu) calculated using the tight-binding model 
assuming that the bandgap is zero, (e): the four bands of bi- 
layer graphene in the absence (left) and in the presence (right) 
of the bandgap, with the interband transitions shown with ar- 
rows, (f): Re G{uj) calculated assuming that the bandgap Ag 
is present as given by the red solid curve (Ref.[T7|). 

The assignment of the optical conductance structures 
to interband transitions is given in FiglJI^d). Apart from 
the discussed strong peak structures B and C there is 
an onset-like structure A which corresponds to a transi- 
tion between the low-energy bands 2 and 3, which has 
the same origin as the onset-like structure observed in 
monolayer graphene 0]. The onset frequency is twice 
the Fermi level with respect to the Dirac energy, which 
is in bilayer graphene proportional to \Vg — Vgo\ with 
a coefficient determined by 70. In the measured spec- 
tra (Fig[3Jb)) we observe such a structure showing the 
same (within the experimental uncertainty) dependence 
on the gate voltage. This confirms that 70 is close the 
value used in the calculation (3.12 eV). This observation 
is in accordance with a recent measurement of Li et al. 
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[241 ■ Interestingly, there is a second onset- like structure, 
with the onset energy showing a similar V-shape depen- 
dence on the gate voltage but shifted with respect to the 
structure A by about 71. The structure is due to the 
onset of transition D (1 ^ 3) for the electron doping 
and transition E (2 ^ 4) for the hole doping. There is a 
significant enhancement of Re G{uj) close to the 'vertex' 
point Lu « 71, V ~ Vgo where the two onsets are close to 
each other 0, Q ■ One can clearly see a similar structure 
on the experimental graph. Thus the tight binding model 
reproduces most of the features of experimental spectra. 

Now we address the issue of the gate-induce bandgap 
A„ between the low-energy electron and hole bands 
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Its manifestation in the infrared spectra 
was first calculated (assuming that 73, 74 and A = 0) in 
Ref.i. In FigHf) we show the result of a calculation 
where we keep the all aforementioned band parameters 
and add a gate-dependent difference in electrostatic po- 
tential between the two planes. We use a curve Ag{Vg) 
from Reffirl. shown as a red line in Fig[3l^f), where the 
charge screening effects were treated self-consistently. We 
assume, as it was also done in Reflrzl. that contaminant 
molecules shifting the charge neutrality point away from 
Vg = act as an effective top gate electrode. In this case 
the bandgap vanishes not at Vg = Vgo but at Vg = ~KjO- 
At the highest gate voltages of our experiment the gap 
value is expected to be of the order of 0.1 eV. 

According to the calculation, the opening of the 
bandgap indeed brings some extra features to the spec- 
tra. All of them are due to the flattening of bands 2 and 
3, as shown in FiglSlJe), which results in a strong increase 
of the density of states of these bands. The first feature 
(marked A') is an enhancement of the optical intensity of 
the transition 2^3. Although this enhancement largely 
shows up at photon energies below the experimentally ac- 
cessible region, its tail spreads up to about 0.2 eV. The 
second feature is the appearance of high-frequency satel- 
lites (marked E' and D') to the peak- like structures B 
and C. These satellites correspond from transitions 2^4 
and 1 — > 3 respectively. The energy separation between 
the central frequencies of peaks B and E' as well as be- 
tween C and D' is close to the energy of the bandgap 
and could be therefore read directly from the conduc- 
tance curves. Note that the interband structures A', E' 
and D' involve the same band pairs as the structures A, 
E and D respectively. However the former ones are exclu- 
sively due to transitions within a very small momentum 
region around the Dirac point. 

We notice that experimental spectra (Figl3]^b)) show 
an enhancement of conductance similar to the high- 
frequency tail of the structure A'. However the satel- 
lite structure E' and D' are not obviously present in the 
data. Thus the tight binding model that is quite suc- 
cessful in describing the main infrared features, is only 
in partial agreement with the data as far as the bandgap 
related features are concerned. This fact is perhaps the 



largest surprise of our study. We can only speculate 
about the possible reasons. First of all, the satellite fea- 
tures might be smeared out by doping inhomogeneity, 
due to the flake corrugation, contaminant molecules or 
other factors. However, the calculation already takes a 
large broadening (about 0.02 eV) into account. A second 
possibility is that the actual bandgap is smaller than the 
prediction of a simple model that does not take into ac- 
count interaction effects, so that the satellites E' and D' 
cannot be easily separated from the main peaks. A third 
possibility is that the gap can be partially filled with im- 
purity states (2^. Future experimental and theoretical 
developments are required to resolve this intriguing is- 
sue. 
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